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Abstract. 

I report on the current status of QCD phase diagram at vanishing baryon density. I focus on 
the QCD phase diagram with three degenerate quark flavor using Highly Improved Staggered 
Quarks on A^r = 6 lattices. No evidence of a first order phase transition in the pion mass 
window of 80 < < 230 MeV is found. The pion mass at the critical point where the chiral 
first order phase transition ends is estimated to be m'^ < 45 MeV. 



1. Introduction 

Due to the asymptotic freedom feature of QCD, it is expected that the QCD matter will undergo 
a phase transition from hadronic phase to the quark gluon plasma phase at sufficient high energy 
density. The nature of such a phase transition crucially depends on the values of quark masses 
(light mi and strange m^), the number of quark flavors {Nj) and the baryon density {^ib) [I1[2]EI- 
The conjectured QCD phase diagram at vanishing baryon density from effective theories on the 
basis of the renormalized group together with the universality [3] is summarized in the left plot 
of Fig. [TJ As sketched in Fig. [H in the pure gauge case, i.e. in the limit of ^ oo and — )■ oo, 
QCD has an exact Z(3) symmetry and the phase transition is of first order. In the chiral limit 
of Nf = 3 case, i.e. rai = rus — >• 0, the phase transition is also known to be a first order. In the 
chiral limit of a Nf = 2 theory, i.e. mi ^ and m<j — )• oo, if Ua_{^) symmetry is restored the 
relevant symmetry is 0(2) x 0(4) and the transition is likely first ordeiH; if Ua_{^) symmetry 
is broken the relevant symmetry becomes isomorphic to that of the 3-d 0(4) spin model and 
the transition becomes second order belonging to that universality class. In the intermediate 
quark mass region, there is no real phase transition but a (rapid) crossover takes place from the 
hadronic phase to the quark gluon plasma phase. All the first order phase transition regions 
are separated from the crossover by lines of second order phase transitions belonging to the 3-d 
Ising Z{2) universality class. Thus the first order region for Nf = 3 case, the second order 0(4) 
line for Nf = 2 case and the second order Z{2) line are supposed to meet at a tri-critical point 
characterized by a certain value m''" of the strange quark mass. For the system with massless 
Nf > 3 fiavors, the order of its chiral phase transition is always first order and is independent 
on the fate of Ua{1) symmetry. 



^ We discuss implicitly the case with number of colors TVc = 3 in this proceedings. 

^ A second order phase transition is also allowed with a different symmetry breaking pattern [3j. 



The QCD phase diagram at vanishing baryon density has been studied extensively on the 
lattice. It has been well established that QCD with Nf = has a first order phase transition 
from a confined phase to a deconfined phase at the transition temperature Tc ~ 270 MeV [5]. 
With physical values of the quark masses m„ = = m^^^^ and rrig = mFj^^^ it is found from 
lattice simulations with staggered fermions that QCD does not have a true transition but a rapid 
crossover and the pseudo critical transition temperature Tpc ~ 154 MeV [6]. The location of the 
physical point with respect to the tri-critical point is not fully determined. More specifically, it is 
not yet clear whether mf > m*^* or mFj^^^ = m}" or mPj^^'^ < m*^* . If mpj^^'^ > m}" then in the 
limit of ^ one should see a second order transition belonging to the 3-d 0(4) universality 
class, if mFs^^^ = m^P the tri-critical point is a Gaussian fixed point of the 3-dimensional (j)^ 
model and its critical exponents take the mean field values if mFj^^^ < m*''* then in 

— 7- limit first one should cross through a second order transition belonging to the 3-d 
Z{2) universality class and then end up in the first order transition region. Recent lattice QCD 
calculations with improved staggered fermions indicate that QCD transition in the chiral limit 
of mi with rUg = mFj^^'^ is likely of second order and belongs to 0(4) (0(2)) universality class [8j. 
To confirm this scenario, further study with better discretization scheme of QCD action and the 
investigation on the fate of Ua{^) symmetry at finite temperature [9] are crucially needed. 




Figure 1. Left: schematic QCD phase transition behavior for different choices of quark masses 
(^u,d) "T-s) at zero chemical potential. Right: The critical surface swept by the chiral critical line 
at finite chemical potential. A QCD chiral critical point may exist if the surface bends towards 
the physical point. The left and right plots are taken from Ref. [10] and Ref. [11], respectively. 

For QCD with Nf > 3 a first order phase transition has been observed at small quark mass 
on coarse lattices using unimproved [12j as well as the p4-improved staggered fermions |13]. 
However, the value of pion mass at the critical point where the first order phase transition ends, 
m^, is not yet determined in the continuum limit and its value at finite lattice cutoff varies from 
~300 MeV to ~50 MeV, strongly depending on the lattice spacing and discretization schemes 
of the action [TUHHIITI]. The value of and the curvature of the surface swept from the 
chiral first order phase transition region at finite may be relevant to the (non-)existence of a 
QCD chiral critical point as sketched in the right plot of Fig. [1] although other possibilities for 
generating a second order transition at the physical values of quark masses have been discussed 
in Ref. [15]. Recent lattice QCD studies based on the imaginary chemical potential approach 
disfavor the existence of such a QCD chiral critical point [16j. This result, however, is obtained 
from simulations using unimproved staggered fermions on coarse lattices and may be prominently 
infiuenced by lattice cutoff effects. Thus it would be interesting to investigate this topic using 
an improved fermion action on finer lattices. As a starting point to determine the curvature of 



the chiral critical surface, the value of needs to be determined. In this talk we thus focus on 
the QCD phase diagram with three degenerate Highly Improved Staggered Quarks at vanishing 
baryon density and discuss the value of m^. 



2. Universality class near critical lines 

As mentioned above in the chiral limit of three flavor QCD the transition is of first order. 
This first order phase transition becomes weaker with finite quark mass and ends at a critical 
quark mass rric sitting on a second order phase transition line. The universal properties of this 
critical point belong to a Z(2) symmetry. The proper order parameter could be a mixture of 
two relevant quantities, e.g. a combination of the chiral condensate with the pure gauge action 
Sg [TQ] , M = + i^Sg) oc (m — rric)^/^ and the susceptibility of the order parameter 

M, i.e. xm/T"^ 0^ {m — mcy/^~^ . Here for Z{2) universal class 1/5 — 1=-0.785. To estimate 

the value of rric, we ignore the contribution from Sg and will use the quark chiral condensate as 
an approximate order parameter for the second order phase transition at the critical point. 
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Table 1. Parameters of the numerical simulations. 



As mentioned in Ref. [lOj the chiral first order phase transition region from lattice QCD 
simulations shrinks with reduced cutoff effects. This is also observed by the lattice QCD 
calculation using the stout action jTl]. The Highly Improved Staggered Quark (HISQ) action 
achieves better taste symmetry than the asqtad and p4 actions at a given lattice spacing [T7]. 
Simulations presented here have been carried out with three degenerate quark flavors using the 
HISQ action. The quark masses vary from 0.0075 to 0.0009375 corresponding to pion masses in 
the region of 80 < ^ 230 MeV. In this proceeding we update the results reported in Ref. |10) 
by performing additional simulations with am = 0.00375 on 16^ x 6, 12^ x 6 and 10^ x 6 lattices 
and simulations with quark mass am = 0.025 on 24^ X 6 lattices. 



3. Results 

We show the volume dependence of quark chiral condensate with am = 0.0009375 (m.,^ ~ 80 
MeV) in the left plot of Fig. [21 The volume dependence at high temperature is small while 
at low temperature it is relatively large. This is obvious since at high temperature the scale 
of the system is controlled by the temperature T while at low temperature the scales are the 
chiral symmetry breaking scales. In the right plot of Fig. [2] we show the volume dependence of 
chiral susceptibility. If the system with quark mass am = 0.0009375 is in the first order phase 
transition region, the chiral susceptibility should scale with volume. However, such volume 
scaling behavior is not observed in the system with pion mass mj^ ~ 80 MeV. 




Figure 2. Volume dependences of chiral condensates (left) and chiral susceptibilities (right) 
with quark mass am = 0.0009375. 




Figure 3. Left: time history of the quark chiral condensate near /3c with am = 0.0009375 on 
24^ X 6 lattices. Right: Chiral condensates as a function of /?. 



In the left plot of Fig. [3l we show the time history of quark chiral condensate near the 
pseudo critical (3 value for our lowest quark mass on 24^ X 6 lattices. There is no evidence for 
the coexistence of two phases. We then investigate the temperature dependence of the chiral 
condensate at different quark masses shown in the right plot of Fig. [3l No evidence for a 
discontinuity in {'4>^l^) as function of /3 at all quark masses is found. Together with the evidence 
from Fig. [2] we conclude that there is no first order phase transition in the system with pion 
mass ranging from 230 to 80 MeV. 

In the left plot of Fig. H] we show the disconnected part of the chiral susceptibility as a 
function of /3 at all quark masses. The pseudo critical temperature decreases with decreasing 
quark mass. It is expected as the hadronic degrees of freedom in the system become lighter and 
thus they become more easily excited in the thermal heat bath. They then can contribute to 
the energy density of the system and thus trigger the onset of a phase transition already at a 
lower temperature. The peak height of chiral susceptibility grows with decreasing quark mass as 
the system is approaching the first order phase transition region. We then performed a scaling 
fit using an ansatz of A{m — rric)^"^/^ to the inverse heights of chiral susceptibility peaks. The 
results are shown in the right plot of Fig. HI The intercept of the solid line with the x axis gives 
an estimate for the critical quark mass. Taking into account that in the infinite volume limit the 
inverse peak height becomes larger (see Fig. [2] left) an upper bound for the value of the critical 
pion mass at the critical point can be obtained, i.e. m^ < 45 MeV. 

4. Conclusion 

We have performed 3-flavor QCD simulations using the HISQ action on Nj- = 6 lattices with six 
pion masses in the region of 80 < mj^ < 230 MeV. Through the study of quark chiral condensates 
and chiral susceptibilities, we found no evidence of a first order chiral phase transition in this 




P ma 

Figure 4. Left: the disconnected part of chiral susceptibility as a function of /3. Right: the 
scahng fit to the inverse height of chiral susceptibility peaks. The filled points denote the inverse 
peak heights of chiral susceptibilities on largest volume at all quark masses while the open points 
are rescaled from the filled points by the volume effects. The solid line and dashed line are fitting 
curves using Z(2) critical exponents. 

pion mass window. The upper bound on the pion mass at the critical point is estimated to be 
around 45 MeV by Z(2) universal scaling analysis considering the finite volume effects. In the 
quark mass plane this upper bound for the critical point in the 3-flavor QCD is thus at about 
(mpj^y/270, mpj^y/270), which is far away from the physical point at (mpj^y/27, m-pj^y). Together 
with the results from Ref. [H], our results suggest that the first order phase transition region 
is very small and thus the critical surface swept by the chiral critical line at finite chemical 
potential has to bend towards the physical point with a very large curvature to affect the nature 
of the transition at physical values of quark masses at a small chemical potential. 
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